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ABSTRACT 



We report the discovery of gamma-ray pulsations from the nearby isolated millisecond pulsar 
PSR J0030-I-0451 with the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope 
(formerly GLAST). This discovery makes PSR J0030-t-0451 the second millisecond pulsar to be 
detected in gamma-rays after PSR J0218-I-4232, observed by the EGRET instrument on the 
Compton Gamma Ray Observatory. The spin-down power E = 3.5 x 10'^^ ergs s~^ is an order of 
magnitude lower than the empirical lower bound of previously known gamma-ray pulsars. The 
emission profile is characterized by two narrow peaks, respectively 0.07 ± 0.01 and 0.08 ± 0.02 
wide, separated by 0.44 ± 0.02 in phase. The first gamma-ray peak falls 0.15 ± 0.01 after the 
main radio peak. The pulse shape is similar to that of the "normal" gamma-ray pulsars. An 
exponentially cut-off power-law fit of the emission spectrum leads to an integral photon flux above 
100 MeV of (6.76 ± 1.05 ± 1.35) xlQ-^ cm'^ s'^ with cut-off energy (1.7 ± 0.4 ± 0.5) GeV. 
Based on its parallax distance of (300 ±90) pc, we obtain a gamma-ray efficiency L^/E 2± 15% for 
the conversion of spin-down energy rate into gamma-ray radiation, assuming isotropic emission. 

Subject headings: pulsars: general - pulsars:individual (PSR J00304-0451) - gamma-rays:observations 



1. Introduction 

Two distinct pulsar populations are the "nor- 
mal" and "millisecond" pulsars (MSPs), the latter 
being rapidly rotating neutron stars (P < 30 ms) 
with very small period increases {P < 10~^^ s/s). 
MSPs represent roughly 10% of the pulsars listed 



in th e ATNF online catalogue (jManchester et al 



20051 ) . In the classical framework of the magnetic 



braking model, MSPs are old stars with charac- 
teristic spin-down ages r — P/{2P) ~ (0.1 - 10) 
X 10^ yrs, and characteristic surface dipole mag- 
netic fields, B.urf ^ 3.2 x lO^^ (PP)^/^ < 10^° 
G. Most MSPs are in binary systems. They are 
thought to have been spun-up by the accretion of 
matter and thus transfer of angular momen tum 
from a binary companion (jAlpar et al.l 119921 ). In 
some binary systems, the companion is evaporated 
by the strong rela tivistic wind produced by the 
millisecond pulsar ( Ruder man et al. 1989f ). which 
becomes isolated, like PSR J0030-F0451. 

PSR J0030-f0451 was discovered by two in- 
depend ent radio sur veys, the Arecibo Drift Scan 
Search (ISomeill200n[) and the Bolog na sub millisec- 
ond pulsar survey lD'Amicoll2000h . Its spin-down 
age T is 7.6 x 10^ yrs. The analysis of radio tim- 
ing residuals showed a significant annual parallax 
of 3.3 ± 0.9 mas, lea ding to a distance m easure- 
ment of 300 ± 90 pc ( Lommen et al. 2006[ ). useful 



for lu minosity estimates. The Cordes fc Laziol 
(l2002l) model of Galactic electron distribution 
predicts a distance of 317 pc, in good agree- 
ment with_thB_2a£allax_measurement. In addi- 
tion, iLommen et al. I (l2006l) also argued that the 
Shklovskii effect (jShklovskilllSTOf) on the pulsar's 
first period derivative P is less than 1%, hence 
rejecting the possible contamination of proper 
motion in P and reinforcing the determination 
of the spin-down energy rate E. This P — 4.87 
ms and P = 1.0 x 10^^" s s~^ pulsar hence has 
an E = 47r2/(P/p3) of 3.5 x 10^^ ergs s-\ tak- 
ing the moment of inertia / to be lO''^ g cm^. 
The magnetic field strength at the stellar surface 
Bs = 3.2 X 10^^ G y/pP for this pulsar is 2.2 x 
10® G, and its characteristic age r = P/{2P) is 
7.7 X 10^ yrs. 

This pulsar was detected in X-rays with ROSAT, 
during the final days of the mission, sho r tly af - 



ter its discovery in radio (Becker et al. EooO). 



The MSP was then observ ed by XMM-Newton 
([Becker fc Aschenbachll2002r) . Both telescopes re- 
vealed a broad X-ray pulsation profile, with a 
pulsed fraction compatible with 50% and the two 
X-ray peaks separated by 0.5 in phase. Unfortu- 
nately neither of the two X-ray telescopes could 
provide the X-ray alignment relative to the radio 
at the time of the observations, due to the lack of 
accurate clock calib ration. We address thi s align - 

■^Slnstitucio Catalana de Recerca i Estudis Avangats (ICREA), Bdai3elJj6iEi3S\JSiri)elow. iBecker fc AschenbachI (|2002l ) 
■^^Consorzio Interuniversitario per la Fisica Spaziale (CIFS), I-101Sll(3teeidLd,hla41;^he X-ray Spectr um is Cons i stent w ith 
^^Dipartimento di Fisica, Universita di Roma "Tor Vergata", I-OObMligJipaa'fcll^iythermal. iBogdanOV et al. ( 2008( ) in- 
**School of Pure and Applied Natural Sciences, University of Kalmar, SE-391 82 Kalmar, Sweden 
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voked the presence of a hydrogen atmosphere to 
model the X-ray spectrum. 

Six pulsars were detected in gamma-ra ys with 



high confidence with the EGRET telescope (jThompson et 
19991) . and more recently AGILE and LAT re- 
ported the detect ion of pulsed gamma-rays from 



PSRJ 2021-h3651 (iHalnern et al.ll2008l : lAbdo et al 

The LAT also discovered a young pul- 



sar in the supernova remnant CTAl (jAbdo et al 
20081 ). All eight are normal pulsars. Gamma- 



ray pulsations with 4.9 a statistical significance 
were reported for th e 2.3 ms pulsar J0218-I-4232 



(|Kuiper et al.l l2004f ) with EGRET. No pulsed 



gamma-ray em ission had previous ly been detected 



from an MSP (jPierro et al.lll995l ). LAT observa- 



tions confirm the detecti on of gamma-ray e mission 
from PSR J0218-h4232 (|Abdo et al.ll2009cl ). 

The revise d EGRET catalogue of ga mma- 
ray sources of ICasandjian fc Grenieii (|2008l ) lists 
the new source EGR J0028-I-0457 with no 3EG 
counterpart, with a 95% confidence contour ra- 
dius of 0.51°. The millisecond pulsar PSR 
J0030-h0451, for which iHarding et all (|2005l ) pre- 
dicted a gamma-ray flux superior to that of the 
marginal EGRET detection PSR J0218-f4232, is 
located 0.5° from EGR J0028+0457 and was hence 
suggested as a possible counterpart. The launch of 
the Fermi observatory, formerly GLAST, on 2008 
June 11 provided a new opportunity to study the 
emission of gamma-rays by MSPs. 

This paper describes the discovery of pulsations 
from the isolated pulsar PSR J0030-h0451 in the 
Fermi Large Area Telescope (LAT) data, making 
this the first high-confidence detection of a mil- 
lisecond pulsar in gamma-rays. We also discuss 
the timing analysis of yjsHM-Newton data with ac- 
curate clock calibration, providing the phase align- 
ment of the radio. X-ray and gamma-ray emission 
relative to each other. 

2. Observations and analysis 



Th e LAT instrument is described in lAtwood et al 
( 20091 ). Gamma-rays convert to electron-positron 
pairs in the tracker, consisting of tungsten foils 
interleaved with layers of silicon microstrip detec- 
tors. The tracker provides direction information. 
Below the tracker is the calorimeter in which par- 
ticles interact in cesium iodide crystals, giving 
most of the energy information. The detector is 



surrounded by the anticoincidence detector, which 
helps reject the charged cosmic-ray background. 
The LAT is sensitive to photons with energies be- 
2§v 20 MeV to over 300 GeV. Its large field of view 
of 2.4 sr, large effective area of 8000 cm^ on-axis at 
1 GeV, improved angular resolution (0.5° of 68% 
PSF containment at 1 GeV for events collected 
in the 'front' section with thin radiator foils), 
the scanning observing mode, and small trigger 
deadtime of 26.5 /xs, make the LAT much more 
sensitive than EGRET. Ground tests using cosmic 
ray muons demonstrated that the LAT measures 
event times with a precision better than 1 /is. On 
orbit satellite telemetry indicates comparable ac- 
curacy. The software timing chain from the GPS- 
based satellite clocks through the barycentering 
and phase-folding software has been shown to be 
accur ate to better than a few /zs ( Smith et al.l 



20081) 



PSR J0030-t-0451 is a very stable pulsar. The 
timing solution used for this object has been de- 
rived from observations made with the Nangay 
radiotelescope, near Orleans, France. About 
seven hundred observations starting in July 1999 
were used, contemporaneous with the Fermi LAT 
dataset and bracketing the XMM observation. Be- 
tween 1999 and 2002, data were recorded with the 
Navy Berkeley Pulsar Processor (NBPP) back- 
end. This instrumentation was designed and built 
at the Naval Research Laboratory in coUabora- 
tion with the University of California , Berkeley 
(|Backer et al.lll997t lobster et al.lll996[) . The sys- 
tem covers 1.5 MHz per channel for a total band- 
width of 144 MHz centered at 1360 MHz. In 
timing mode, the data were folded for 15 min- 
utes over 256 bins covering the full period of 
the pulsar PSR J0030-f0451 (with a resolution 
of ~ 20 /Lts). Since late 2004, TOAs are ob- 
tained from the coherent pulsar instrumentation 
currently in use at Nanqay, the Berkeley- O r leans- 



Nanqay (BON) backend (jTheureau et all l2005t 



ICognard &: Theurcau 20061 ) ■ The coherent dedis- 
persion is performed within 4 MHz channels over 
a total bandwidth of 64 MHz (128 MHz since July 
2008) centered at 1398 MHz. Compared to those 
recent high quality data, the old NBPP TOAs 
were degraded due to the absence of dedispersion 
inside the 1.5 MHz channels (also producing large 
systematics when the pulsar scintillates). The 
bulk of radio observations were done at 1.4 GHz 
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(1360 ± 72 MHz before 2002, 1398 ± 32 MHz 
after 2004 and 1398 ± 64 MHz after 2008). 2 
GHz observations were made to better constrain 
the dispersion measure, necessary for the ahgn- 
ment of pulsar profiles at different wavelengths. 
The TEMPCQ package was used to build a timing 
solution from the recorded radio times of arrival, 
determine d through a standard cross-correlation 
procedure (|Tavloilll992l ). The mean time of arrival 
uncertainty is 3.6 /xs for the recent BON observa- 
tions (the bulk of the data), while this is 8.6 fis for 
all the dataset. The post-fit rms is 3.7 /its, fitting 
for the pulsar's frequency and first-derivative, and 
taking its proper motion into account. The tim- 
ing parameters used for phase-folding gamma-ray 
events are given in Table [1] the timing residuals 
after the fitting procedure are plotted in Figure 
[TJ Doubling the error in the parallax given by 
TEMPO, and adding in quadrature an additional 
0.3 mas un certainty due to the s olar wind, as ex- 
plained in iLommen et al. we measure a 



parallax of 4.1 ± 0.7 mas, which is consistent with 
the Lommen et al value. The dispersion mea- 
sure we derive is also consistent with the value 
of 4.3 328(3) pc cm~^ quoted in iLommen et al 
(l2006l) . 



The LAT data considered here were taken dur- 
ing Fermi's first-year all-sky survey, starting 2008 
August 3, through November 2. The pulsar is 
located well outside of the Galactic plane (1 = 
113.141° and b = -57.611°) and is hence in a re- 
gion of low Galactic background. Aiming to have 
a good signal-to-noise ratio over a broad energy 



^http: //www. atnf . cslro . an. research/pulsar/tempo | 




Fig. 1. — Timing residuals as a function of time 
for the model given in Table[TJ Data between 1999 
and 2002 were recorded using the NBPP backend, 
the BON backend was used after 2004. 



range, we used an energy-dependent region of in- 
terest of 6* = Max[0.9 - 2.1 Logio(EGov), 1-5] de- 
grees around the pulsar position. A larger fraction 
of the PSF is included at high energies, where 
there is little background contamination. Since 
Fermi was operating in survey mode, the contribu- 
tion of the Earth's gamma-ray albedo to the back- 
ground was negligible. Finally, the "diffuse" class 
events were kep t. A description of e vent classes 
can be found in Atwood et al. I (l2009l) . After ap- 
plication of the above cuts, we obtained a dataset 
of 563 events over 100 MeV. 

PSR J0030+0451 was observed by XMM- 
Newton on 2001 June 19-20. The observations 
with the pn camera spanned 29 ks, but a soft pro- 
ton flare affected approximately 8.8 ks of the ex- 
posure. The data were reduced using Version 8.0 
of the XMM-Newton Science Analysis Softwar^. 
This version solves the timing problem. The pn 
camera was used in timing mode with the thin fil- 
ter. The data were reduced using "epproc" . The 
event lists were filtered, so that - 4 of the pre- 
defined patterns (single and double events) were 
retained, as these have the best energy calibration. 
We used the data between 0.3 and 2.5 keV as this 
had the best signal-to-noise. The event times were 
converted to Barycentric Dynamical Time, using 
the task "barycen" and the same coordinates as 
used for the Fermi data. The XMM-Newton abso- 
lute timing accuracy is found to be as good as 300 
fi^ or ~ 0.06 rotations of PSR J0030-f0451. In 
order to obtain such fine temporal resolution with 
the pn in timing mode, the CCD is read out conti- 
nously, causing the events for the target source to 
be smeared out in the Y-direction. To extract all 
the events from the source, we used the standard 
procedure of creating a one dimensional image, 
by binning all of the raw data in the Y-direction 
into a single bin. The spectrum was extracted 
using a range of 7 pixels centered on the pulsar, 
in the X-direction. The background spectrum was 
extracted from a similar neighbouring region, free 
from X-ray sources. 



^ http : //xmin . vilsp a . esa . es/sas/ 1 

• ^http: //xmm2 . esac ■ esa. lnt/docs/docunients/CAL-TN-0045-1-0 .pdf | 
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Table 1 

Timing ephemeris for pulsar PSR J0030+0451. 



Parameter Value 



Right ascension, a 00:30:27.4303(5) 

Declination, 5 04:51:39.74(2) 

Proper motion in right ascension, (mas yr^"^) -5.3(9) 

Proper motion in declination, fi§ (mas yr~^) .... -2(2) 

Parallax (mas) 4.1(3) 

Epoch of position determination (MJD) 52079 

Pulse frequency, 1/(3"^) 205.530699274922(9) 

First derivative of pulse frequency, u (s^^) -4.2976(4) X 10^^° 

Epoch of ephemeris (MJD) 50984.4 

MJD range 51343 - 54757 

Number of TOAs 651 

Rms timing residual {^s) 3.69 

Dispersion measure, DA'I (cm~^pc) 4.333(1) 

Solar system ephemeris model DE200 



Note. — Figures in parentheses are the nominal Icr TEMPO uncertainties in 
the least-significant digits quoted. Epochs arc given in TDB units. 



3. Results 

3.1. Gamma- ray and X-ray profile analysis 

Figure [2] shows the phase histogram of the 
events with energies greater than 100 MeV in 
panel a, along with the reference radio profile used 
to derive the timing ephemeris in panel c. The 
value for the phase histogram shown in panel a 
of Figure [2] is 121 for 29 degrees of freedom, in- 
dicating that the probability that the pulsation 
is actually a statistical fluctu ation is 3.1 x 10~^^ . 
The bin independent 77-Test (|de Jager et al.lll989h 
gives a value of 123. The derived chance oc- 
curence probability of the null hypothesis, non- 
pulsed emission, is below 4 x 10^^. The null hy- 
pothesis is hence ruled out, for a single trial. 

The zero of phase in the reference radio pro- 
file is defined to be at the maximum of the first 
Fourier harmonic of the signal, transferred back to 
the time domain. The maximum of the main radio 
peak's second subpulse is at 0.036 in phase. The 
gamma-ray pulse profile comprises two peaks (see 
Figure [21 panel a). There is a shift between the 
first gamma-ray peak (PI), occuring at 0.15 ± 0.01 
in phase, and the main radio peak. The gamma- 
ray peaks of PSR J0030+0451 are very sharp. We 
fit PI and P2 with two-sided lorentzians, to take 
into account the different widths for the leading 
and trailing edges. For the first peak, the fit 
gives a full width at half maximum of 0.07 ± 0.01 



in phase, that is, 340 /is. The second gamma- 
ray peak, P2, lags PI by A$ = 0.44 ± 0.02 in 
phase. The fit of the structure between 0.45 and 
0.6 places the peak at 0.59 ± 0.01 in phase with 
a fuU width at half maximum (FWHM) of 0.08 ± 
0.02 in phase. 

LAT phase histograms in two energy bands are 
given in Figure [3l with 30 bins per rotation. The 
pulsar is faint in the 100 to 500 MeV band. In this 
energy band, P2 is prominent relative to PI: tak- 
ing PI (resp. P2) between 0.1 and 0.25 in phase 
(resp. 0.45 and 0.60), the P1/P2 ratio is found 
to be 0.47 ± 0.17. The upper panel of Figure [3] 
shows the events with energies over 500 MeV. In 
this band PI dominates P2, with a ratio of 1.64 ± 
0.48. There hence seems to be a spectral depen- 
dence of the gamma-ray profile. Within error bars, 
the P1/P2 ratio seems to increase as a. function of 



energy , conversely to the Vela pulsar (jAbdo et al 



2009al ). The first gamma-ray peak hence seems 
harder than the second one. More photons are 
needed to perform phase-resolved spectroscopy. 
We note that in the 100 MeV to 500 MeV energy 
band, P2 seems broader and closer to PI than it 
does for energies greater than 500 MeV. However, 
we cannot be conclusive about this trend due to 
the low number of photons in P2: more statistics 
might reveal bridge emission between PI and P2, 
in which case P2 would be narrower than the 0.08 
+/- 0.02 value quoted above. We also note that 
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Fig. 2. — Multi- wavelength phase histograms 
of PSR J0030+0451. Two pulsar rotations are 
shown, a: Gamma-ray phase histogram of PSR 
J0030+0451 at E > 100 MeV, within an energy- 
dependent ROI. Each bin is 0.033 in phase, or 160 
IIS. b: 20 bin 0.3-2.5 keV XMM-Newton phase 
histogram, c: Radio profile obtained at 1.4 GHz 
used to build the ephemeris. 



there is no evidence for pulsed emission from 20 
MeV to 100 MeV. Extrapolation of the hard, ob- 
served spectrum (see spectral analysis below) pre- 
dicts only a few photons in this energy band, and 
any pulsation is almost certainly obscured by the 
rapidly rising background below 200 MeV. 

The overall gamma-ray emission profile is 
reminiscent of the younger pulsars Vela, Crab, 
Geminga, PSR B195H-32 or PSR J2021+3651, 
especially over 500 MeV. The phase separation of 
~ 0.4 between the two peaks is a common fea- 
ture ( Thompson 20041 ). The pulse profile also 
shows a main gamma-ray peak lagging the main 
radio component by 0.15 ± 0.01, similar to the 
main gamma peak lagging 0.11 to 0.16 after the 
radio pulse for Vela, PSR B195H-32 and PSR 
J2021-h3651. The radio to PI and PI to P2 sepa- 
ration are in agre ement with the outer-gap puls e 



profile model of iRomani &: Yadigarog 



which compares the gamma-ray peak separation 
A with the radio to gamma-ray lag 6 for a given 
magnetic inclination angle a. This is also in agree- 
ment with the two-pole caustic and slot gap mod- 
els for gamma-ra y emission (Wattcrs ct al. 200^ 



Dvks et al.ll2004l). With radio p olarization mea- 
surements, Tomme^^r^ (2000) found the mag- 



netic inclination angle a to be most probably 62°, 
that is a mostly orthogonal configuration. 

Panel b in Figure [2] shows the 0.3 to 2.5 
keV XMM-Newton phase histogram for PSR 
J0030-I-0 451. The ove r all pr ofile is similar to that 
found bv lBecker" et al.l (|2000^ and' Becker fc AschenbachI 



(l2002[ ). 1 -hough their observations did not yield the 
X-ray to radio alignment because of inaccurate ab- 
solute timing. Here the ~ 300 /is absolute timing 
accuracy represents ± 1 bin in panel b. The X-ray 
and radio are hence consistent with being phase 
aligned. This result supports the idea that the X- 
ray and radio emission have common origins in the 
pulsar magnetosphere, and that gamma-rays are 
produced in a different region. In that sense, PSR 
J0030-f 0451 is different f r om P SR J0218-f4232: 
according to iKuiper et al. ( 2000l ). the 2.3 ms pul- 
sar has its radio. X-ray and gamma-ray compo- 
nents aligned. 
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3.2. Gamma-ray spectral analysis 




Fig. 3. — 30 bin phase histograms for PSR 
J0030+0451 with the LAT, in different energy 
bands. Two rotations are shown. The horizontal 
dashed Une shows the background level estimated 
from a surrounding annulus. 



Using the standard maximum-likelihood spec- 
tral estimator gtlike in the Fermi Science Toolfl 
we performed a spectral analysis of the gamma 
sample. The diffuse Galactic background and ex- 
tragalactic emissions are taken into account, as 
well as the instrument response, which is a func- 
tion of the photon energy and incidence angle rel- 
ative to the telescope axis. The dataset used for 
the spectral analysis spans the same time inter- 
val as for the timing analysis, but this time we 
retain events within 15° of the millisecond pulsar. 
Because of present uncertainty in the instrument 
response below 200 MeV, events with energy be- 
low 200 MeV are rejected. 

Figure [4] shows the phase-averaged differential 
energy spectrum. The corresponding power-law 
with exponential cutoff modeling the data is given 
by: 



'dE 



E 



103 MeV 



,-E/E^ 



(1) 



In this expression, E is in MeV, the prefac- 
tor term iVo = (1.84 ± 0.38 ± 0.37) xlO"" 
ph cm~^ s^^ MeV~^, the power-law index F = 
(1.4 ± 0.2 ± 0.2) and the cut-off energy Ec = 
(1.7 ± 0.4 ± 0.5) GeV. The first error is statis- 
tical, the second is systematic, dominated by dif- 
ferences between the a priori expectations of the 
Fermi LAT effective area and the on-orbit instru- 
ment response. The magnitude of this effective 
area uncertainty is < 10% near 1 GeV, 20% below 
0.1 GeV and 30% over 10 GeV. Analysis improve- 
ments are underway. 

The stability of the results was tested by fit- 
ting the same dataset with a binned maximum 
likelihood estimator, 'ptlike', which computes the 
photon counts in a point source weighted aperture 
in excess of background counts. The fit results, 
shown in Figure[5]for each energy band, are consis- 
tent with those obtained with 'gtlike' within error 
bars. We have also tried a simple power-law fit to 
the data, of the form dN^/dE = N^{E /IGeY)-'^ . 
The fit using the exponential cutoff functional 
form is better constrained, with a difference in the 
log likelihoods of 10.76. A x'^ interpretation of this 



" ^http : //f erml ■ gsf c . nasa ■ gov/ssc/data/analysls/SAE.overvlew . html | 



value leads to a probability of incorrectly rejecting 
the power-law hypothesis of 3.5 x 10~®. 

Integrating Equation ^ for energies > 100 
MeV yields an integral flux fyiooMeV — (6.76 



± 1.05 ± 1.35) xlO-« p h cm-^ s 



The re- 



vised EGRET catalogue of ICasandiian fc Grenier 



(120081) quotes a flux of (10.4 ± 3.1) x 10"® ph cm'^ s" 



for EGR J0028-f0457, based on the summed 
EGRET dataset. Both fluxes are statistically in 
agreement. The pulsar is in a low EGRET ex- 
posure region, mainly detected at large off-axis 
angles where systematic uncertainties in the effec- 
tive area were large. 

The energy flux is Fobs = (4.91 ± 0.45 ± 0.98) 
xlO"" ergs cm"2 s^^ over 100 MeV. The lu- 
minosity of a pulsar can be written as = 
4:7TfnFobsD^ , where D is the pulsar distance, /ji is 
a correction factor containing information about 
the beaming geometry, a.nd Fobs is the observed 



phase-averaged energy flux. Watters et al. ( 20081 ) 
have computed pulse profiles and fn corrections 
for young pulsars. The outer magnetospheres for 
MSPs should be scaled-down analo gues of those of 
young pulsars. If we use the maps in lWatters et al.l 
(|2008l ) for a pulse separation A = 0.4 and a mag- 
netic inclination a — 62°, as inferred from the ra- 
dio polarization data, we find a reasonable match 
for an old, high efficiency two-pole caustic/slot gap 
model for a viewing angle C ~ 75°. In turn this 
implies /n ~ 0.8. This pulse width is not natu- 
ral in a high efficiency, large gap width outer gap 
model, suggesting that the lower altitude two-pole 
caustic picture is a better match to the data, un- 
less the true pulsar efficiency is < 0.1. Adopting 
fn = I for PSR J0030-I-0451, we obtain an ef- 
ficiency rj — L^/E = 15% for the conversion of 
spin-down energy into gamma-ray emission. 

Arond (|l996l ) noted that for the EGRET pul- 



sars, r] is, in good approximation, inversely pro- 
portional to the open field line voltage = 4 x 
-[^Q20p-3/2pi/2^ proportional to V ^, and also pro - 
portional to the open field current ( Hardinj|l981 ). 
With an open field line voltage of 1.2 x 10"'^^ volts 
and an efficiency of 15%, PSR J0030-h0451 seems 
to break from the trend, which would predict a 
higher value. However, the efficiency law may sat- 
urate at lower spin-down rates. Other low E pul- 
sar detections with Fermi should help constrain 
how these old "recycled" pulsars convert their en- 
ergy loss rate into gamma-ray luminosity. In ad- 



dition, the discovery of pulsed emission from PSR 
J0030-f0451, which has a smaller spin down rate 
than previously known gamma-ray pulsars, low- 
ers the empirical E threshold for gamma-ray emis- 
sion by an order of magnitude. This suggests that 
many low E pulsars might be detectable by the 
LAT. 

Attempts to describe the high-energy emission 
from MSPs were recently made based on the two 
main classes of theoretical models: the polar cap 
(PC) and the outer gap (OG ) models. In the PC 
description bv lHarding et al. I lf2005l . charged par- 
ticles are accelerated along the open field lines 
near the magnetic poles to high altitudes. The 
high-energy spectrum consists of three main com- 
ponents. The emission of photons up to 100 MeV 
is dominated by synchrotron radiation from elec- 
trons. Over 100 GeV, photons are produced by 
inverse Compton radiation from electrons. Curva- 
ture radiation from electrons dominates the pho- 
ton spectru m between 1 and 100 GeV. For PSR 
J0030-H0451. lHarding et~all (|2005l) provide a pre- 
diction of 4.25 GeV for the curvature radiation 
cutoff energy. This value differs from the 1.7 GeV 
we find in this analysis. In their model, the ex- 
pected curvature radiation flux over 100 MeV is 
Fcr{> lOOMeV) - 3 X lO^^ph cm'^ s'^ for 
our pulsar. Our measured integral flux disagrees 
with t heir expectation. However, iHarding et al 



(|2005[ ) computed sky-averaged spectra; more de- 
tailed three-dimensional sums may be needed to 
make predictions for an individual viewing an- 
gle. On the other hand, in the OG descrip- 
tion of the high - energy radiation from MSPs by 
Zhang fc Chengl (2003), a strong multipole mag- 
netic field exists near the stellar surface. X-rays 
are then produced by the backflow current of the 
outer gap, which is a vacuum gap close to the light 
cylinder. These X-rays consist of a non-thermal 
power-law component, plus two thermal compo- 
nents. Gamr na-rays are produced in the outer 
gap. Though IZhang fc Ch eng (2002) provide no 
application of their model to the case of PSR 
J0030-I-0451, they predict that a gamma-ray MSP 
would be an X-ray MSP. This is indeed the case 
for both PSR J0030-h0451 and PSR J0218-h4232. 
Secondly, they predict that if the X-ray spectrum 
is dominated by thermal emission, which is indeed 
the case for PSR J0 030-I-0451. consistent with be - 
ing purely thermal ( Becker fc AschenbachI 2002h . 
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the gamma-ray emission can extend to ^ GeV 
gamma-rays only, as we observe. 

4. Conclusion 

We described the detection of the millisecond 
pulsar PSR J0030+0451 in gamma-rays using the 
Fermi LAT. Wc now have high confidence that 
there are gamma-ray emitters among millisecond 
pulsars. This provides a new tool for studying 
the magnetospheres of energetic pulsars. De- 
tection of more millisecond pulsars in the LAT 
data may invite revisiting the possible contri- 
bution of unresolved millisecond pulsars to the 
overall Galactic diffuse emission in t he ga mma- 
ray band: for instance IWang et al.l (|2005l ) pro- 
posed that the Galactic centre might contain a 
few thousand unresolved MSPs, contributing to 
the diffuse spectrum detected by EGRET, which 
shows a break at a few GeV. One might also 
expect cumulative gamma-ray emission of mil- 
lisecond pulsars in globular clusters, such as 47 
Tuc, which is thought to contain up to 60 MSPs 
(|Camilo fc Rasiol l2005l ) . The EGRET search for 
emission from globu lar clusters only provided 



1994: Fierroetal 



uppe r limits (,Michelson et al 
19951 ). The Fermi LAT offers new opportunities 



to search for emission from globular clusters in 
gamma-rays. 
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Fig. 4. — Phase-averaged spectral energy distribu- 
tion for PSR J0030+04.51, along with the best-fit 
power law with exponential cut-off. Statistical er- 
rors are shown. 
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